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We derive a microscopic transport theory of multiterminal hybrid structures in which 
a superconductor is connected to several spin-polarized electrodes. We discuss the non- 
perturbative physics of extended contacts, and show that it can be well represented by 
averaging out the phase of the electronic wave function. The maximal conductance of a two- 
channel contact is proportional to (e^//i)(ao/-D)^ exp [— where D is the distance 
between the contacts, ap the lattice spacing, ^{lo) is the superconducting coherence length, 
and Lo* is the cross-over frequency between a perturbative regime (w < lo*) and a non per- 
turbative regime {lo* < w < A). The intercontact Andreev reflection and elastic cotunneling 
conductances are not equal if the electronic phases take a fixed value. However, these two 
quantities do coincide if one can average out the electronic phase. The equality between the 
Andreev and cotunneling conductances is also valid in the presence of at least one extended 
contact in which the phases take deterministic values. 



I. INTRODUCTION 



Transport of correlated pairs of electrons in multiterminal configurations has recently focused an important in- 
terest. One possible line of research is motivated by the possibility of creating entangled pairs of electrons from a 
superconductor This may lead to fundamental tests of quantum mechanics in solid state, or to new ways of 

manipulating quantum information. In a different context, the interplay between superconductivity and magnetism 
offers novel functionalities in multiterminal devices : the various transport channels occurring at several neighboring 
superconducting- ferromagnet interfaces depend in a subtle way on spin polarizations and geometry p|-p^. At the 
theoretical level, there is a need for a transport theory in multiterminal hybrid structures involving superconducting 
and spin-polarized elements. These structures have so far not been the object of any experiment. As a consequence, 
one of the objective of the present time models is to predict what should be measured in future experiments. 

In ferromagnet ~ superconductor junctions, it is well established that Andreev reflection is suppressed by an increase 
of spin polarization Q . This is because Andreev reflection can take place only in the channels having both a spin- 
up and a spin-down Fermi point. This theoretical prediction has been probed experimentally by two independent 
groups [@,||. Spin polarized Andreev reflection 0| and related effects can even be used to measure the Fermi 
surface spin polarization. 

We are concerned here with more sophisticated systems in which a superconductor is connected to several electrodes, 
which can be ferromagnetic or normal metals. In the case of ferromagnetic electrodes, it will be crucial to take into 
account the existence of a very small coherence length. We neglect any diffusive effect in spite of the fact that they 
lead to a rich physics [p^-p^. As a consequence, our models should apply to point contacts having a dimension 
much smaller than the diffusive mean free path. The fabrication of such contacts in multiterminal configurations 
may seem difficult in view of the present day technology, and this is why there are no available experiments on 
these systems. However, there are interesting phenomena taking place in these multiterminal systems. For instance, 
Andreev reflection can become non local. Namely a spin-up electron from a given electrode A can be Andreev reflected 
as a hole in a different electrode B. This effect has been studied theoretically by Byers and Flatte in Ref. for 
normal metals, and in Ref. |^ in the ferromagnetic case which contains the richest physics. It is in our opinion crucial 
to develop the most general theoretical description of this phenomenon. In this respect, two approaches have been 
developed recently. One is based on the analysis of the lowest order processes appearing in perturbation theory 
Another approach is non perturbative, but relies on effective Green's functions B. It is a very natural task to work 
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out the microscopic theory of transport in baUistic multiterminal hybrid structures. Transport theory will be solved 
exactly by means of Green's function techniques [ pT| . We incorporate in our theoretical description two realistic 
constraints: 

(i) Multichannel effects which are expected to play a central role in quantum point contacts involving ferromagnetic 
metals. The radius of the contact can be smaller or larger than the phase coherence length of the ferromagnetic 
metal. 

(ii) The strength of the tunnel amplitude is small in low transparency contacts, and large in high transparency 
contacts. Our approach is non perturbative. Therefore, the tunnel matrix element can take arbitrary values. 
As a result, we can derive transport in the presence of arbitrary bias voltages. 

An ingredient that is not incorporated at the present stage in the model is the reduction of the superconducting gap 
associated to the proximity effect. 

We will use this non perturbative approach to address the following physical questions: 

(i) It has been already established that the multiterminal hybrid system should be described by a conductance 
matrix. The matrix elements encode all information about the current flowing in a given electrode, in response 
to a voltage applied in another electrode. How does the crossed conductance behave when the voltages are close 
to the superconducting gap ? What is the maximal value of the crossed conductance ? 

(ii) The superconductor Green's functions contain not only an information about non local processes, but contain 
also an information about the phase of electron propagation. In extended contacts, there are many phases 
coming into account. Can these phases be considered as random quantities ? 

(iii) There are two propagators associated to a superconductor: the ordinary and the anomalous propagators. One 
can easily realize that after phase averaging, the ordinary propagator is identical to the anomalous propagator 



(see section [II A 3). As a consequence, in the tunnel approach and for unpolarized contacts, the averaged 
Andreev reflection conductance is equal to the the averaged elastic cotunneling conductance Q . Is this identity 
still valid in the presence of large interface transparencies ? 

The article is organized as follows. Some preliminaries are given in section ||. The form of the Green's functions 



is derived in section III. The solution of the model with two single-channel electrodes is presented in section IV 
The general solution with an arbitrary number of single-channel electrodes is presented in section As a particular 



are solved in section VII 



example, we discuss in section VI the physics of a model with three single-channel electrodes. Multichannel electrodes 



II. PRELIMINARIES 



A. The Keldysh method 



We will use Green's functions techniques to solve transport theory. There is an advanced (G^), retarded (G^) and 
Keldysh (G"*"'") Green's functions p8| , p^ . Each of these Green's functions is a 2 x 2 matrix in Nambu representation. 
The Dyson equation for the advanced and retarded Green's functions takes the form 



The Dyson equation for the Keldysh component is given by 



G+- = 



I + G 



(1) 



(2) 



Eqs. 1^, H are written in a compact notation in which the convolution involves a summation over time variables and 
space labels. S is the self energy, which contains all couplings present in the tunnel Hamiltonian. The notation g is 
used for the Green's functions of the disconnected system (i.e. with S = 0) while G refers to the Green's functions of 
the connected system (i.e. with S 7^ 0). We will use the notation 



{t,t') = 



fA,R^t,t') g^-«(t,t') 



(3) 
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for the Nambu representation of the advanced and retarded Green's functions, with 



We win also denote hy p — ilni(g"^) the Nambu representation of the density of states: 



P 



Pa Pf 
Pf Pa 



(4) 
(5) 

(6) 



with pg = ^l'ni{g^) and pf = ■i-Im(/"^). The Nambu representation of a given timncl matrix element connecting sites 
a and a is ta.aO"^, where is one of the Pauli matrices. We use a notation in which the "sites" of the superconductor 
are represented by the Greek symbols a, (3 , 7... The sites in the normal metal electrodes are represented by the Latin 
symbols a, b, c , ... The explicit form of the Keldysh Green's function connecting the two sides of a given interface is, 
from (0) 



(7) 



The strategy is first to use (|^) to calculate the advanced and retarded Green's functions and next use (0) to calculate 
the Keldysh Green's function. The current can be obtained easily from the Keldysh Green's function |M : 



dio 



(8) 



The spin- up (spin-down) current is given by the 11 (22) matrix element of the Nambu representation. 



B. A useful trick on the spectral current 

The systems of interest here are made of a single superconductor connected to an arbitrary number of external 
normal metal electrodes. It turns out that there exists some tricks that can be used to simplify the calculation of the 
current in such systems. One of these tricks is the following. 

The terms in the first summation in Eq. ^contain a prefactor nF(w — /is) because Ui belongs to the superconductor. 
The terms in the second summation contain a prefactor nF{t-o — ^Mai) because ai belongs to a ferromagnetic electrode. 
The trick consists in realizing that the terms containing nF{uj — ps) coincide exactly with the terms containing 
npioj ~ fJ-oi)- As a result, the total current can be written as an integral over energy of the spectral current: 

4 = ^ du) [nF{uj - Pat) ~ npiuj - ps)]1k,i{^)- (9) 



III. FORM OF THE GREEN'S FUNCTIONS 



In this section, we present a derivation of the form of the Green's functions that will be used throughout the 
remainder of the article. This will give us the opportunity to discuss the relevant parameters of the model. 



A. Green's functions in the superconductor 

1. Effective Green's functions 

In some cases, it will be useful to describe the superconducting and ferromagnetic reservoirs in terms of effective 
Green's functions. It was already shown by one of us in Ref. Q that effective Green's functions can be used to construct 
a consistent non perturbative "toy model" version of transport theory. In this approach, the superconductor is viewed 
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as zero dimensional: its dimensions are shorter than the coherence length. The superconducting effective Green's 
function takes the form IIitI: 



-{OJ- ^s)2 



(10) 



and we will consider in the following the limit 77 — > 0. The Keldysh component is given by g'^^ [uj) = 2iTrnp{u! 
/is)p(w), with p{uj) — -lni{g^) the density of states. The ferromagnetic electrodes are described in a similar way: 



Pis 

P2,2 



(11) 



where and p2,2 are respectively the spin-up and spin-down densities of states. 



2. Spectral representation 



To address more realistic models, it is useful to restore the dependence of the Green's functions upon space co- 
ordinates. The Green's function are evaluated in terms of a spectral representation. Let us start with the pairing 
Hamiltonian 



H = £kci cr + Aid ci + AfeC.- tCr 



(12) 



with k = \k\. We use the notation ■fg = eg — Mi with — h'^k'^/{2m) for the kinetic energy. After standard 
manipulations, the spectral representation is found to be standing for ai,aj): 



UJ - {ps + Ek) + IT] UJ - {ps - Ek) + ii] 
1 1 

UkVk 



UJ - {fis + Ek) +iri UJ - {ps - Ek) + IT] 



where M the number of sites in the superconductor. We used the standard notation 

Ek - VA2 + {ikf 
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(13) 
(14) 

(15) 
(16) 

(17) 



for the quasiparticle energy, and the electron and hole coherence factors. 



3. Form of the Green's function 



The spectral representation (13), (|T^ is valid in any dimension. We now restrict our discussion to the case of a 
three dimensional superconductor. We first perform the angular integration and next use the residue theorem to make 
the radial integration. Note that it is crucial to carry out first the angular integration. This ensures the existence of 
well defined convergence properties when we use the residue theorem to make the radial integration. The final result 
is 



9i,3 



1 



h? 2'K\Xi ■ 



■ exp 



sm (fi 



VA2 -{uj- fis f 



2e(c.) 



(18) 



-(cj-//s)±i77 
A 



A 




'10' 




{uJ - Ms) ± _ 


— cos tp 


1 


1 
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with = kF\xi — Xj \ and ap is the length of the elementary cell. The coherence length appearing in ( |T^ ) is 

1^ +00 it > A. 



(19) 



for the zero- frequency coherence length, with ep the Fermi energy. In the case of 



We used the notation ^(0) 

It — 7I R A 

two point contacts a and b treated explicitly in section [V, the Green's function provides a coherent coupling 
between charge transport at the two contacts. We end-up this section with three remarks. First, we note that with 
cosiyj = the Green's functions arc identical to the effective Green's functions given in section IIIAl. Second, we 
recover the usual free-fermion Green's function in the limit w ^ A: 



1 



,1 



ti^ 2T:\xj 



Finally, we will discuss in detail the role played by phase averaging. Using the notation ((...)) = / one can show 
that (((ffi.j)^)) = (((/i,j)^))- This identity implies that in the tunnel limit the average Andreev reflection conductance 
is equal to the average elastic cotunneling conductance. 

B. Green's functions in the ferromagnetic electrodes 

The Green's function in the ferromagnetic electrodes are diagonal in Nambu space. The form of the Green's function 
is taken as 



91,0 A"^) = -'^ 



1 



?i2 2tt\ 



X 1 X n 



■ exp 



exp 



I "^"i "^3 I 



21 



(20) 



where the phase is given by — k^p^ \xi — Xj\ 
the spin- up and spin-down Fermi wave vectors: 



and Z^'^'' is the phase coherence length. There is a mismatch between 



k 



'2m 



and /lox is the exchange field. At some point, it will be convenient to assume that the phase takes the particular 
value Lp^'^^ = 0. With this special value of the phase, the form of the 3D Green's function ( pO| ) is identical to the 
effective Green's function (^l|). The coherence length in ferromagnetic metals is much shorter than in usual metals 



so that ferromagnetism is often treated in a semi-classical description (see [20-22 
Aharonov-Bohm oscillations reported in Ref 



For instance, the absence of 



]12| shows that the coherence length in Co is smaller than 0.3/im. This 
can be incorporated in our model by considering that the "dissipation" 77 is not a small parameter. This results in a 
finite coherence length, which is spin-dependent, and inverse proportional to the strength of dissipation: 



ri v2to 



This simple phenomenological model contains the relevant physics associated to ferromagnetic metals. For instance, 
the phase coherence length of spin-up electrons is larger than the spin-down coherence length. The Green's function 
( pO| ) is infinite when Xi = xj, which is also the case for the superconductor Green's function (p^. Local quantities 
can be obtained by using \xi — Xj\ — ao instead of Xi = Xj. With this condition, the local density of states of the 
ferromagnet is given by 



P - 2^2 21 



Spin- up electrons have thus a larger density of states than spin-down electrons. 



IV. SINGLE CHANNEL ELECTRODES: (I) TWO ELECTRODES WITH 100% SPIN POLARIZATION 



In this section, we consider a model in which a single channel spin-up electrode a and a single channel spin-down 
electrode h are in contact with a superconductor. We assume in sections IV A and IV B that a is a spin-up channel 



and 6 is a spin-down channel. The result for parallel spin orientations is given in section IV C 
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a spin-up electrode 



b spin-down electrode 



FIG. 1. Schematic representation of the model considered in section [V Two ferromagnetic electrodes are in contact with a 
superconductor. A voltage Vs is applied on the superconductor while the voltages Va and Vt are applied on the ferromagnetic 
electrodes. 



A. Derivation of the transport formula 

1. Solution of the Dyson equation 

Let us first calculate the Nambu representation of the propagators. The starting point is the chain of Dyson 
equations given by (||): 

/^a.a "I r ^a.a "I f T^a.a T^a.h ~\ f /^a.a 

(21) 



- Qa,a ■ 
















Qb,a 







+ 




j^b,b 







where we used the notation K^'^'^^ = -ai,anai,ai-ai,a,-j-a,,a, ^ rpj^g solution of Eq. M is 



Ga,a y^a^h 

1,1 "-^1,2 

^b,a y^b,b 

•-^2,1 '-^2,2 



1 



where the determinant I?af is given by 



1 - 



1 - \t''^\ht29''A - r^''?\t''^?9Z9t2r'^'^''^ 



(22) 



(23) 



and g and / have been defined as the components of the Nambu matrix in (^ and (|l^). If not specified, all Green's 
functions in a given formula stand as well for advanced and retarded functions, and similarly for the determinant 
I?AF given by Eq. ( |2^ ) and the determinant 2?f that will be introduced latter. The matrix in Eq. contains the non 
vanishing Nambu components of the renormalized propagator. Because we assume a complete spin polarization, the 
other Nambu components are vanishing. For instance G°^\ — C?2'i = G'^\ = 0. 



B. Exact expression of the current 

Using the expression of the Keldysh propagator (see Appendix ^), we deduce the final expression of the spin-up 
current in electrode a: 



a, a „a,a 

a 



It'i = -47r^|<Q,ap J duj [npiuJ - /-ia) - npic^ - /is)] pI'.iP 

1 - 52,2 9 l-|4,/3|ff2,2 9 



+ 2iT7\ta^a\'^\tbJ3f J did [n^ (w - fla) - (w ^ Ps)] pT,i92 



a, a _b,b.A 
2 



1 



ja,fS.A j(3,a,A 



1 u \2 b,b,R l3,fJ.R 
1 - Kb,/3| 52.2 9 



a.a b,b.R 



2iTT\ta,a\ \tb,f3\ duj [npiLo - fia) - npiuj - ps)] Pil92 



(24) 



(25) 



(26) 
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-2|, |2 



u \2b,b,A f3,[3,A 
Kb,/3| 92,2 9 



(27) 



which generahzes the result obtained by Cuevas et al. in Ref. in the case of a single conduction channel. From the 
density of state prefactors, we see that there are two type of contributions to the current: (i) The quasiparticle current, 
which is proportional to the product of the density of state in the superconductor (pg) and one of the ferromagnetic 
electrodes (for instance pa); (ii) The crossed Andreev current which is proportional to the product paPb of the density 
of state in the two ferromagnetic electrodes. The term ( |24|) contributes only to quasiparticle current. The term ( [27| ) 
contributes only to Andreev reflection. The mixed terms (p5|) - (|2^) contribute both to the quasiparticle and Andreev 
current. 



C. Two-terminal conductance matrix 



To understand the meaning of the transport formula ( 
structure in terms of a differential conductance matrix: 



|) , it is useful to describe transport across the multiterminal 



where the matrix elements are given by 



Qa ,a Qa,b 
Gb.a Gb.b 



din 



Ga.,a,{Va,Vb) = ^(K,Vfa) 



(28) 



(29) 



The conductance matrix (^8|) encodes all information about transport in the two-terminal structure. The off-diagonal 
matrix elements should satisfy a symmetry relation: Ga,biVa, Vb) = Gb.aiVb, Va)- If the electrodes have an antiparallel 
spin orientation, subgap current is transported by Cooper pairs if < A, in which case we have la ^ lb- This implies 
an additional symmetry relation: Ga,a{Va,Vb) = Gb,a{Va,Vb), and Gb,b{Va,Vb) = GaA^a,Vb). If the electrodes have 
a parallel spin orientation, subgap current is due to elastic cotunneling, in which case /„ = — /{,. The additional 
symmetry relation reads Ga,a{Va,Vb) = -Gb,a{Va,Vb), and GbA^^'^b) = -Ga,b{^a,Vb). 



1. Sub-gap conductance matrix: effective Green's functions 



In this section as well as in section [VC2, we assume that costp = so that we can use effective Green's functions 
to evaluate the transport formula and work out the basic physics on simple grounds. The validity of this assumption 



will be discussed in section [V D 



Superconductor 



/ / S 
/ / 
/ / 



FIG. 2. The diagram associated to Andreev reflection in the two-channel model. 



a. Antiparallel magnetizations: The sub-gap current of the two-channel model with antiparallel magnetizations 
originates from the non local Andreev reflections which are shown schematically on Fig. |. This can be seen by 
inserting the effective Green's functions into the transport formula: 
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^^'^=^^'^=- \v.^%\^ [/-^(^=Hr, (31) 

where = 7r|fa,aPpa is the spectral line-width associated to electrode a, and a similar expression holds for Fj. We 
used the fact that g*^'^ and /"■'^ are real numbers below the superconducting gap. The expression of 2?af is the 
following: 

|2?AF(a;)|2 = {1 - F'^F^ - (/"./3)2] f + (r-5".« + r^/./3)' . (32) 

h. Parallel magnetizations: The same calculation can be done if the electrodes have a parallel spin orientation. 
We find 

= = - \vZ%? = """^^ ' ^''^ 

with 

|©f(w)|^ = {1 - FT* _ (5«>/')2] }^ + (r»5«'« + F''/''^)^ . (35) 

There are two differences between the situations with antiparallel and parallel spin orientations. First, the Andreev 
reflection transport with antiparallel spin orientations is controlled by the anomalous propagator /"■'^ while the 
elastic cotunneling transport with parallel spin orientations is controlled by the ordinary propagator The second 
difference is in the sign of the off-diagonal conductance matrix elements. The four matrix elements have the same sign 
in the case of Andreev reflection because transport is mediated by Cooper pairs. The off-diagonal matrix elements 
have a sign opposite to the diagonal matrix elements in the case of elastic cotunneling because transport is due to 
single electron tunneling between the two electrodes. 



2. Conductance matrix above the superconducting gap: effective Green's functions 

Let us now assume that the voltage Va is above the superconducting gap and that the clcetrocies have an antiparallel 
spin orientation. With the notation g^-P'^^^ = ±i|(7"'^|, and /"'^'^-^ = ±i|/"''^|, the extra diagonal terms of the 
conductance matrix take the form 



yb,a — ^2 \J I ' 



which should be evaluated at the energy w = 14. The diagonal conductance matrix element is the sum of a crossed 
contribution and a quasiparticle contribution: Qa^a = —Gb,a — Ga^a- The quasiparticle contribution is the sum of a 
direct and a crossed term: 



Gl^a = (1 + r.l/'^l)^ - rl\r'^\' , (36) 



, 2 



where the denominator T>af is a real number: 

Vaf = = pfp = 1 + r-^r" - ir-^H + ri<7"'"| + r'lg^'^l. 

In the limit |;ro — — * +00 in which the separation between the contacts becomes very large, the extra diagonal 
conductance matrix elements are vanishingly small. The quasiparticle term reduces to the conductance of a single 
channel metal - metal contact: 

^ 47rFaPg 
{l + ra\g"'"\f' 
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D. Phase resolved versus averaged conductance 



Given the form (18) of the superconductor Green's function, we see that the conductance depends exphcitly on the 



electronic phase difference tp — fc^la^a — Xb\- This leads us to calculate the conductance in two different ways: 

(i) The phase-resolved conductance G{f)- We will focus more especially on the case Lp = 7r/2. For this special value 
of the phase difference, the Green's function of the superconductor coincides with the effective Green's function 
in the limit \xa — Xb\ ^ a^. 



(ii) The averaged conductance 



^ /.27r 



((^(^)» - ^ (37) 



This phase averaging is used to mimic the physics of extended contacts that will be considered later in section VI] 



To determine the role played by phase averaging, we compare the phase-dependent and average conductances (see 
the end of section III). It is visible on Fig. ^ that the effective Green's function conductance (i.e. with ip = it/2) 



follows closely the average conductance. This shows that the effective Green's function conductance contains already 
the relevant physics, as far as the two-channel problem is concerned. 




-15 -10 -5 

Log^-ro) 

FIG. 3. Variation of the logarithm of the crossed conductance \ogQa,b versus log(A — lo) {lo < A). Qa,b is in units of e^//i. 
The points correspond to the phase-resolved conductance with ip = tt/2, namely, to the generalized effective Green's functions. 
The solid lines correspond to the average conductance (|37|). We used the parameters m = 0.01, kp = 1, ep = 50, ao = 1, 
A = 1. The distance between the contacts is _D = 100 and the superconductor coherence length is = ef/CcfA) — 50. 



It is also visible on Fig. g that there is a cross-over energy uj* (thus a cross-over voltage V* = (^* /e) at which 
the crossed conductance reaches a maximum, li lo < uj* the crossed conductance behaves like Ga.b ~ 1/(A — w). If 
Lu* < Lu < A, the crossed conductance behaves like Ga,b ^ A — Only when lu < lu* does our approach coincide with 
the lowest order tunnel perturbation theory. The analysis based on generalized effective Green's functions is simplified 
if one assumes that D ^ uq, but D can be small or large compared to the superconducting coherence length ^{lo) 
(see Eq. p^ ). The behavior in the energy range w* < < A is non perturbative, and can be understood by retaining 
in Daf only the leading divergence, which is generated by the quartic terms (see Eq. ^2[): 



|2?af|' 



(/ 



a,p\ 



from what we deduce the expression of the crossed conductance 



r 4 fao^2 (2Titi 



A2 



exp 



D 



(38) 



valid in the energy range u> > u)* . The expression of the cross-over energy to* is obtained by equating the quadratic 
and quartic terms in 2?af- This leads to 
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, ma? 



A,/i-(r3). (39) 



To obtain the conductance in the energy range u) < uj* , we expand Daf up to order 



where g^°'^ — g"'" = g^'^ denotes the local propagator in the superconductor. One can show easily that f"'^ <^ g^°^ 
because D qq, from what we deduce 

,o o / ma?, \ ^ , , 

|PAFp=.l + 2r2 -JI -— — . (40) 



Vaf is close to unity only when tu < u>o, with 



^0 = , =■ (41) 



Comparing ( p9D and (41), we see that uj* is larger than ujq but luq and lu* have the same order of magnitude is T is 
small (which is the case on Fig. |3|). We deduce that I?af = 1 in the relevant energy range uj < luq < lo* . From what 
we obtain the conductance in the energy range lu < loq: 

identical to the one obtained in lowest order perturbation theory. Evaluating (|3^ ) and ( ^2|) at uj — uj* or oj = ujo leads 
to the maximal value of the conductance: 

where the numerical prefactor of order unity cannot be obtained from this simple estimate. 

It is well known from the BTK scattering approach |^ that the conductance of a normal metal - superconductor 
contact is equal to twice the quantum of conductance e^//i per spin channel if a; = A, regardless the value of the 
interface scattering. The same behavior occurs in the Keldysh formalism treatment by Cuevas et al. [|l^ . This type 
of resonance can be properly described only with a non perturbative approach. Eq. |4^ constitutes a generalization 
of the BTK behavior to the case of spatially separated contacts having a phase difference ip = tt/2. Given that the 
average and phase-resolved conductances follow closely each other (see Fig. ^), it is expected that ( ^ ) is also valid 
for the average conductance, but with an extra reduction factor. 



V. SINGLE CHANNEL ELECTRODES: (II) GENERAL SOLUTION 

Now we consider a model in which N spin-up electrodes and M spin-down electrodes are in contact with a su- 
perconductor (see Fig. ^). All the necessary details can be found in Appendix The final form of the current 
is 

N 



fc=2 •' 



2m^{-f+^t''^^°'H^>"'^>' dtj[nF(w-/iai) -?^F(t^-M5)] X (46) 

k=2 •'l^"! 
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a I spin-up electrode 
• 

spin-up electrode 



b| spin-down electrode 



hyi spin-down electrode 



FIG. 4. Schematic representation of the model with N ferromagnetic spin-up electrodes and M ferromagnetic spin-down 
electrodes. 



1 ( rai,t3k,R CaR CaA _ fai.,l3k-A f^A CaR \ 
N 

7 O 



k=2 

ai,ai afc.Ofc -■ai,Qfc,i?,~Qfc,ai,A 
M 



fe=l 



(47) 



(48) 



ai,ai 6fe,fcfc fai,l3k,R fPk.ai.A 

where we used the notation 



f 



a., ft _ 



(49) 



(50) 



for the renormahzed propagators. There are three types of processes involved in the transport formula: (i) The 
quasiparticle term (ji^) which is proportional to p"^'"^pg; (ii) The elastic cotunneling term (^^ in which spin- 
up electrons from electrode ak are transfered into electrode 1. The elastic cotunneling terms are proportional to 
P±A^^ pV'i^'^ (iii) The Andreev reflection term (48) which are proportional to pW""^ P22'' ■ The mixed terms ( ^5|) and 
( P^ ) contribute the three types of processes. 



VI. SINGLE CHANNEL ELECTRODES: (III) THREE ELECTRODES WITH 100% SPIN POLARIZATION 

Let us now consider a three-terminal problem. We consider that each of the three electrodes has 100% spin 
polarization (see Fig. ^). The aim is to have Andreev and cotunneling processes occurring in the same multiterminal 
device, which allows a direct comparison of these two basic processes. The transport formula can be deduced easily 
from the general solution obtained in section M. 



A. Three-terminal conductance matrix 



The three-terminal conductance matrix generalizing ( Pq ) takes the form 



Q = 



Oa.a Oa^b 
Qb,a Gb,b 
Gc,a Gc,b 



Ga.c 
Gb,c 
Gc,c 



(51) 



Let us assume that electrode a has a spin-up magnetization, and electrodes b and c have a spin-down magnetization 
(see Fig. ^). We use the same procedure as in section EVC 1 to obtain the conductance matrix elements. Namely, we 
assume that cos = and replace the Green's functions by effective Green's functions. 
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a spin-up electrode 



b spin-down electrode 



c spin-down electrode 



FIG. 5. Schematic representation of the three channel model considered in section VI 



1. Conductance matrix below the superconducting gap: effective Green's functions 



Let us first give the form of the off diagonal matrix elements: 



Ga.b — Gb,a — — 
Qa.c Qc.a 
Gb.c = Gc,b = 



4rjb 

X>Ax>R 



p2 ^f...Pfan _ g^,agP.n^ 



(52) 

(53) 
(54) 



From the signs of these matrix elements, and from the type of propagator involved, we see that Ga,b and Ga,c correspond 
to Andreev reflection while Gh^c corresponds to elastic cotunneling. 




FIG. 6. Variation of the logarithm of the Andreev reflection conductance Qa,c and elastic cotunneling conductance Gt.c versus 
reduced energy cu/A. The points correspond to Andreev reflection and the solid lines correspond to elastic cotunneling. We 
have assumed that site a coincides with site (3 (see Fig. |^. We used effective Green's functions with (p = tt/2. The distance 
between the contacts is D = 100. On purpose, we did not show the behavior for oj ~ A (see Fig. ^). The parameters are 
identical as on Fig. ^ 



2. Conductance matrix above the superconducting gap: effective Green's functions 

Using the notation f^'^ — ±«|/|, and g^'^ — ±i\g\, we obtain the conductance matrix elements above the super- 
conducting gap: 

Ga.b = Gb,a = [\r^\ + Te (!/"' V'^l " I/"' V'^l)] ' (55) 

Ga,c = Gc.a = [ir--^! + n (irV'^'i - ir-V-^i)]' (56) 

Gb,c = Gc,b = [l/'l + (I/' V'"l - l/^'^T'l)] ' . (57) 
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The diagonal coefficient Qa.a takes the form 

ga,a - -Ga.b - ^a,c + [-npgM^'^-^^M^^-^'' 

+ niir^Pf [IZ-'^I + r, (2 + r,|.g^'^|) (|/^V'''I - \G^^^\')] } 

..... ... 

and similar expressions can be obtained for Qi, i, and Gcc- 



2rfcr,|r''^r'V'^l [1 - nvc (l/V'^l 



'^l-l/'^l')]} 



(58) 

(59) 
(60) 
(61) 



B. Phase resolved versus average conductance 



Let us now consider the special situation in which a ferromagnetic electrode with 100% spin polarization is at a 
distance D away from a normal metal electrode having a zero spin polarization. Namely, we assume that sites a and 
/3 coincide (see Fig. p[). T his provides the simplest model containing both Andreev reflection and elastic cotunneling. 

Following section IV D, we evaluate the phase-resolved conductance with ipa,b — fa.c — fb,c = tt/2. With this 



particular value of the phases, the phase-resolved conductance coincides with the effective Green's function conduc- 
tance. We have shown on Fig. I the energy dependence of the Andreev reflection and elastic tunneling conductances, 
evaluated with ip = 7r/2. It is visible that the Andreev reflection conductance is larger than the elastic cotunneling 
conductance by a factor A/lu. The Andreev conductance coincides with the elastic cotunneling conductance only 
when uj is close to A (see Fig. This behavior can be understood from the effective Green's function conductance 
(Eqs.H-H). 




log(A-a)) 

FIG. 7. The same as Fig. n but for energies li; ~ A. 



To obtain the average conductance, we assume that the phase variables are independent random variables, and 
average the conductance over all possible values of the phases: 



m) = 



27r 



27r 



27r 



■Gifa.b, fa,c, Vb,c)- 



(62) 



It is visible on Fig. || that after phase averaging, the Andreev reflection conductance coincides exactly with the elastic 
cotunneling conductance below the superconducting gap. For energies close to A, we find again the existence of a 
maximum in the conductance at an energy w*. The predictions of perturbation theory are valid in the energy range 
Lo < oj* . In this energy range, all conductance spectra on Fig. ^ can be deduced from each other by a simple rescaling 
(see the insert of Fig. 0) . 



VII. MULTICHANNEL ELECTRODES 
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b 

o 



-12.5 



-17.5 




FIG. 8. Variation of the logarithm of the Andreev reflection (points) and elastic cotunneling (solid lines) conductances 
versus energy uj. We averaged the conductance over all possible phase configurations (see Eq. |62| ). The insert shows the energy 
dependence of the rescaled conductance Q{lu)/Q{lu = 0). The parameters are identical as on Fig. ^ 



o 

3 

-o 
a 
o 
o 

(U 
3 



spin-up 
channels 



spin-down 
channels 



FIG. 9. Schematic representation of the tight binding model. The contacts between the spin-up ferromagnet and the 
superconductor are noted ai,..., a^- The contacts between the spin-down ferromagnet and the superconductor are noted 
bN- 



A. Transport formula 

We want to determine whether extended contacts have a physics identical to the single channel contacts considered 
in sections 0, and 0. We are thus lead to investigate the following situations: 



(i) The phase averaged conductance of extended contacts. There are N{N — l)/2 phases (fij associated with a 
contact having N channels. The phase-averaged conductance is obtained by averaging the conductance over all 
possible values of these phases: 

(ii) The phase-resolved conductance of extended contacts. This is the conductance of an extended contact where 
the phases ipij are deterministic and take the particular value (pij — kp\x.i — Xj| as given by Eq. |l^. 

Let us consider a model in which a multichannel fully polarized spin-up electrode is in contact with a superconductor. 
At a distance I?, there is another fully polarized spin-down electrode. The tight binding model is represented on Fig. |^. 
There is one block "a" made of N fully polarized spin-up channels and another block "b" made of M fully polarized 
spin-down channels. The only difference with section ^ is the existence of a propagator gai,aj, 5bi,6j with i ^ j. 



The form of the Dyson matrix is still given by (Bl). Compared to (B2), there is an additional summation in the 
coefficients of the Dyson matrix: 
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k 
k 

j^ai,aj _ "y gai^ak-j-ak.ak gak.aj-j-aj ,a.j 
k 



X 



bi,bi 



(64) 
(65) 
(66) 
(67) 



The derivation of the transport formula is similar to section For instance, the subgap Andreev conductance is the 
sum of all possible Cooper pair transmissions: 



pG-F„ geFi 



(68) 



where "p G Fa" ("g S Ft,") means that the p runs over all possible channels in electrodes a and b. We used the 
notation 



J/3,, Op _ 



_^ap,bg 



(69) 
(70) 



for the renormalized propagators. The same formalism can be used to handle more complicated situations involving 
an arbitrary number of electrodes having arbitrary spin polarizations. 



2.5 
2 

1.5 
1 

0.5 



Andreev Reflection 
Elastic Cotunneling 





0.25 0.5 0.75 1 

fully 1 down / lup 

polarized unpolarized 

FIG. 10. Variation of the ui — Andreev and elastic cotunneling average conductances as a function of Z'^'/i'^'. The section 
of the two electrodes is circular, with a radius R = 2.1. There are 26 channels in each electrode. The average conductances 
are normalized with respect to the average Andreev conductance with uj = and Z'-^' = 0. The conductance is distributed as 
a function of the phases. The errorbars indicate the root mean square of the conductance distribution. The parameters are 
identical to Fig. ^ 



B. Andreev reflection versus cotunneling 

Let us consider a system in which two multichannel electrodes are in contact with a superconductor: (i) a ferro- 
magnetic electrode; and (ii) a normal metal electrode with no spin polarization. 

Following the discussion in section [II B, we use the ratio Z^^^/Z^^-* to parametrize spin polarization. There is no 
spin polarization if Z^-^^/Z^T) = i and there is a strong spin polarization if Z^^^/Z^T) <^ i \Ye have shown on Fig. |l^ the 
variation of the a; = Andreev and cotunneling phase averaged conductances as a function of spin polarization in the 
ferromagnetic electrode. Elastic cotunneling at w = is vanishingiy small if the ferromagnet is strongly polarized. 
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The average Andreev conductance is equal to the average elastic cotunneling conductance in the absence of spin 
polarization (i. e. with l'^^^ = l^^^). 

Therefore, the phase averaged conductance matrix of the multichannel model behaves like the phase averaged 
conductance matrix of the single channel model. 



C. Extended contacts without phase averaging 



Now we consider extended contacts in which the phases take deterministic values. The phases are given by Eq. 
f = kpl'X-i — We represented on Fig. |ll| the dependence of the Andreev and cotunneling currents as a function 
of ^(^^//(T) for uj = 0, If the number of channels is sufficiently large, we see that the Andreev and elastic cotunneling 
conductances are almost identical for the non magnetic metal {l^^^ /l^^^ = 1). The behavior of extended contacts with 
deterministic phases (see Fig. |l]) is therefore identical to the behavior of contacts with random phases (see Fig. |o|). 
This result, already established in the perturbative regime (tu <C A) is found to be valid at any frequency and barrier 
transparency. We have verified that this behavior is also valid in the non perturbative regime {ut* < oj < A). Notice 
that the same result is found if only one contact is extended, the other one having a few channels. 



t=0.1,lup =1 



I 

§ 
u 




FIG. 11. Variation of the cu = Andreev and elastic cotunneling conductances as a function of i*-^'/?'^' for an extended 
contact with deterministic phases. /'-^'/Z'^' = 1 corresponds to a non magnetic metal. The section of the two electrodes is 
circular, with a radius R — 1.1 (O, 10 channels), R = 2.1 (+, 26 channels), R — 3.1 (□, 58 channels), 7? = 4.1 (x, 98 channels), 
R = 5.1 (A, 178 channels), R = 6.1 (*, 242 channels). The parameters are identical to Fig. ^ The normalizations are identical 
to Fig. 0. 



VIII. CONCLUSIONS 



We have provided in this article a detailed theoretical description of ballistic transport in multiterminal hybrid 
structures involving a superconductor and several spin-polarized electrodes. We have performed a non perturbative 
calculation of the conductance matrix, using Keldysh technique, and focusing on the spin dependence, the geom- 
etry dependence and the energy behavior of the two basic processes : Crossed (intercontact) Andreev and Elastic 
Cotunneling. This generalizes previous perturbative calculations, valid only for low contact transparencies and at 
small voltages. It also generalizes non perturbative calculations made for a single contact and using effective Green's 
functions, which turns out to correspond to a certain choice of the phase (0 = 7r/2) in the electronic propagators. 

A first issue concerns the subgap voltage dependence of the conductances, compared to the usual Andreev conduc- 
tance at a single contact. A maximum is found at a crossover energy (voltage), and the conductance at this maximum 
is reduced compared to the ideal Andreev conductance obtained exactly at w = A in the single contact case. The 
higher the transparency of the barriers, the lower is this crossover energy, below which the perturbative theory is 
essentially valid. 

The other important issue concerns the phase problem. Owing to the different forms of the normal and anomalous 
propagators controlling (in the superconductor) the cotunneling and Andreev processes respectively, the corresponding 
conductances assume different values for single-channel contacts (cotunneling takes place if they have parallel spin 
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polarization, crossed Andreev if they have antiparallel ones) . This was exemphfied here in a three-contact configuration 
with spin polarizations chosen such as both processes can be compared in equivalent geometries. On the other 
hand, averaging the Fermi phase oscillations present in the propagators make the two processes lead to equivalent 
conductance contributions. We have shown that a self-averaging effect occurs when at least one of the contacts has 
many channels. The resulting symmetry (equality of cotunneling and Andreev crossed conductances) was previously 
demonstrated in a perturbative regime, but it now appears as much more general. As a consequence, if at least one 
of the contacts is not spin-polarized, the resulting intercontact conductance is zero, by compensation of cotunneling 
and Andreev processes [Q. Conversely, for spin-polarized contacts this offers a way of measuring, in amplitude 
and direction, the polarization of one contact with respect to the other through the crossed current measurement, as 
proposed in Ref p^ . We believe that multiterminal superconducting-ferromagnet devices have a large, yet unexplored, 
potential in the growing field of spintronics. Further theoretical problems concern the self-consistent calculation of 
the superconducting gap pC| ]. 



APPENDIX A: EXPRESSION OF THE KELDYSH PROPAGATOR: TWO-CHANNEL MODEL 



In this appendix we derive the expression of the Keldysh propagator associated to two single-channel 
with 100% spin polarization (see section IV). We need to calculate the Keldysh component: ta,QG+'~ = 
with i, j £ {a, 6, a, P}, and 

Ka,l3 — ia,a I + C'^^Ja.a 5^'^ il3,bGt,a 

Kp^p = ia.aG'^Jb^pgplph^pGp^ 
K — f c'^ f «+'~ f 4- f r'^ 

Kb,b = ia,aGa,pili.b9h,b tbJjGp^b- 

We also need to calculate ta^aG'^^'^ = J2ij '^i^h 

TC' — C'^ f f ^ i 

^^a.a ~ ^a,a^a,aya,a ^ ' ^a.a'-^a,a 

^'aJJ ^ G^'.a^a.agalp ^/3,fcG'fc^ct 

^'ti,a = G^,bib,pgp^a ^ + ^a,aG'^^a 

^'l3,0 = GaMbJjgp^p ip.bGb^a 

^'b,b ~ Ga\pip,bgt,b tb.pGp ,^. 

The expression of the four terms containing np{ijj — fis) takes the form 



electrodes 



K' — K 



1,1 



A 2 ( \\i |2 O-.O- Oi,OL 

An npl^J - fls)\ta.a\ PlllPg 



1 



-I I, \2 0,0, /I 

1 - \'tb,p\ g2,2 9' 



2b,b,A f3,p,A 



-1 I, \2 0,0, n 



2 b,b,R n,p,R 



K. 



a,p 



A 2 f Ml \2\t \2 a, a a.,P b,b,A pP.a.A 

47r npiuj - iis)\ta,o\ \tb,p\ Pi,iPf 92,2 f ' 

A 2 r Mj. |2u |2 a. a p.a,A b,b.R _pa.P,R 

47r nF[uj ~ p.s)\ta,o,\ \tb,p\ Pi,iPf 92,2 f 



K'p^a - Kp,a 



1,1 



1 



l-lt^'^g'^i-^^-^^'^ 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 
(A6) 

(A7) 
(A8) 
(A9) 
(AlO) 
(All) 
(A12) 

(A13) 
(A14) 
(A15) 
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1 



b,b,A b,b,R pa.f3,R rp,a.A 
i/2,2 92,2 J J ' 



where pg is one of the Nambu components of the superconductor density of states (see Eq. 
/ia and /Xb read 



1 



1 u \2b,b,A B,l3.A\ [1 1, \2 b,b,R I3.I3,R 
1 - Fb,/3| 32,2 9 1 - K&,/3| 52,2 9 



DA-pR 

+ 2iTTnF{uj ~ p.a)\ta,a\'^\tb,0fPlll 
1 



^:^9'2'2^r''^f 1- |4,,P52';^'V-^^^ 



2iTTnF{uJ - ^ia)\ta,a\'^\tb.0fPl]l 



1 



b.b.R fa,0.R r(},a.R 
92.2 J J 



-1 u \2 0,0. J\ 
1 - \'tb,!3\ 92,2 9' 



2 b,b,A B,I3,A 



K'b.b - i'^b,b 



One arrives at the identity 
1 



1 



a,l3,R ff3,a,A 



npiuj - fj.s) 
1 



7 7 ^b,b ~ ^b,b , 

1,1 npy-^ — Mh) L ' J 1,1 



(A16) 

The terms containing 
(A17) 



K,a - Ka,a 



(A18) 

(A19) 
(A20) 



which constitutes for this particular system a proof of the trick on the spectral current (see section 11 B). The 
expression of the Keldysh propagators given in this appendix can be used to obtain the transport formula given in 



section IV B 



APPENDIX B: DERIVATION OF THE TRANSPORT FORMULA WITH AN ARBITRARY NUMBER 

OF SINGLE-CHANNEL ELECTRODES 

We present in this appendix the derivation of the transport formula given by Eqs. (Q) - (|4^), associated to a 
situation where N ferromagnetic spin-up electrodes and M ferromagnetic spin-down electrodes arc in contact with a 
superconductor (see Fig. ||). 



1. Solution of the Dyson equation 

The unknown Green's functions {G"i'"i, G""-"!, G''i'''i, G^^''''^} are the solution of the Dyson equation 



..,G 




,ai 

; ■ ■ 






^ai,ai 




' ga.i,ai 











M 


gfojv.ai 









Qb^,a^ 









QbM.ai 








where the Dyson matrix M. takes the form 



a 



(Bl) 
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where "K"-" is a x TV block, X"-''' is & N x M block The matrix elements of X and Y are 
The solution of the Dyson equation takes the form 



where V is the determinant of the Dyson matrix and Mai,aj are the minors of this matrix. 

2. Solution of the Dyson-Keldysh equation 

To obtain the current through the link ai - ai , we need to evaluate the Keldysh component 



Xai ,Qi /^+.— fai .aiAai .ai.Rfai,ai - + .— 



"ai ,ai /^ai .ai.A 



N 



k=2 

N 



k=l 



Let us start with (]B8|). The first step is to show that 



1,1 



The different terms in this equation are found to be 
and 



fc,Qfc Ql.Qfc KA 

y ■^^iai,ak 



k^l 
M 

k=l 



•^k+N^a^,a^^bk.,0k fai,f3k^^^^^^ 



To evaluate (B9), we first show that 

N 



fc=2 



Using the identities 



(B2) 
(B3) 



^1,1 


(-Y+3 


(B4) 


(-^2,2 


(—Y+j 


(B5) 


'-'1,2 


= — .g'-'^Xb^,.. 


(B6) 


^2.1 




(B7) 



(B8) 
(B9) 

(BIO) 



(Bll) 
(B12) 

(B13) 
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L J 1,1 V / gak,akj) "■^'"■i' 

(_)fe+l 



(B14) 
(B15) 



we obtain 



N 
k=2 



where gai,ak denotes a renormalized propagator. We use a similar calculation to evaluate (BIO) and we deduce the 
transport formula given by Eqs. (H) - (Esf). 
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